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1. Introduction 

It has been known that a large part of sohton equations in 1 + 1 dimensions 
have integrable matrix generahzations, or more generally, integrable multi-field 
generalizations. Such generalizations have been studied for a wide variety of single (or 
two)-component integrable systems by virtue of various approaches such as Jordan 
algebras and Jordan pairs lH, H, H], inverse scattering method 0, |, |, |], Hirota's 
method P, |l^, 0, |T2|, |13|. However, as far as the authors know, matrix generalizations 



of the derivative nonlinear Schrodinger(DNLS)-type systems [0, ^ have not been 
studied throughly. 



In recent papers [[Iq , \L% , Olver and Sokolov made a detailed investigation on the 
DNLS-type systems of the form 

Pt = Pxx + /(-P, Px, Sx) 

St = Sxx + g{Py S, Px, Sx)- 

Here P and S take values in an associative algebra. For simplicity, in the following we 
regard P and 5* as matrix-valued. / and g are non-commutative polynomials of weight 
5, where the weights of dt, dx, P and S are respectively assigned to be 4, 2, 1 and 
1. They made a complete list of the DNLS-type systems (|1.1|) which have one higher 
symmetry of the following form with weight 9, 

Pt Pxxxx ~I~ f i^P^ S, Px: Sx, Pxx J Sxxy Pxxxy Sxxx) 
St Sxxxx ~l~ di^Pj S, Pxi Sxy Pxxi Sxxi Pxxxi Sxxx)' 

Here the commutativity of the two fiows, i.e. dtdrP = drdtP, dtdrS = drdtS, works as a 
strong constraint on the form of / and g in order for the non-commutative polynomials 
/ and g to exist. See ref. |T^, |T^, ^] for a more precise explanation of symmetry 
approach. 

The entries of the list are divided into two kinds. Two systems in the list are 
interpreted as non-abelian analogues of the following integrable system 

Pt = Pxx + 2ap'^Sx + 2appxS - a(3p^s^ 

St = -Sxx + Ws^Px + 2(3sSxP + a(5s^p^ 

for a particular choice of the constants a and (3. The system ( [I.2| ) is linearizable by 
a change of the dependent variables. Thus, we can construct the general solution of 
the system. We often refer to such linearizable systems as 'C-integrable' in Calogero's 
terminology ||2l|. On the other hand, excepting the two entries corresponding to (|1.2|) , 



the scalar-valued counterparts of the systems in the list by Olver and Sokolov are given 
by i 

+ qxx + 4i5gV^ + i(45 - a)gg^r + 5{A5 + a)qh'^ = 
i?"* ~ ^xx + 4i5r qx + i(45 — a)rr^g — 5{A5 + a)r q = 
for special choices of the constants a and 6. This system was generated via a gauge 



transformation for the DNLS equation by Kund |2^ . We can write down a Lax pair for 
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the system ( |1.3| ) with the help of the gauge transformation. According to Calogero's 
terminology, this kind of systems, which are linearizable by the inverse scattering 



formulation, are called 'S-integrable' systems [21 



If a system has one higher symmetry, it has been believed that the system has 
an infinite series of symmetries and thus is completely integrable. However, a system 
proposed by Bakirov was recently proved to be a counter-example to this empirical law 
|23| . Thus, there is no guarantee that the systems in the list by Olver and Sokolov are 
really integrable, although the counter-example seems very exceptional. 

The aim of this paper is to establish the complete integrability of all the matrix- 
valued systems given in [|l^. In a previous paper the authors introduced a Lax pair 



;i-4) 



for the matrix generalization of the Chen-Lee-Liu equation [ p!5|] , 

iQt + Qxx - iQRQ.. = O 
iRt - - iRxQR = O 

which is a member of the list by Olver and Sokolov. In the present paper, we generalize 
the Lax pair for ( |L4|) to be applicable for several matrix systems of the DNLS type in 



The paper consists of the following. In section 2, we extend a notion of the gauge 
transformation to the non-abelian case. We give an explicit expression of the Lax pairs 



for all systems but two in [|17| to show their 'S-integrability'. On the other hand, in 
section 3, we prove the remaining two systems to be linearizable so that they are 'C- 
integrable'. Concluding remarks are given in the last section, section 4. Throughout 
this paper, capital letters denote matrices while small letters represent scalars. Equation 



numbers without 'section-number.' refer to equations in [117 



2. 'S-integrable' systems 

We begin with a brief summary of the matrix generalization of the Chen-Lee-Liu-type 
DNLS equation . A system of linear differential equations for a vector \I' 

= U^ $t = (2.1) 

is consistent if the following equation is satisfied: 

Ut - K- + UV -VU = 0. (2.2) 

The matrices U and V are called Lax matrices or a Lax pair. We introduce the following 
form of the Lax pair. 



u = ie 



-h 



-2/i 



-iR, 



+ c 



Q 



R 



2/2 



+ c 



+ i 
2Q 



O 



iRQ 



:2.3) 



2R 

iQ, + iQRQ 



iRQR 



+ i 



-QR 



RQ 



O 



MRQ.-R.Q) + iRQRQ 



(2.4) 
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Here, ( is the spectral parameter. Ji and I2 are respectively the ni x ni and the n2 x n2 
identity matrices. Q is an ni x n2 matrix and R is an n2 x ni matrix. Putting ( |2.3| ) and 
( p.4|) into (|2.2|) , we obtain the matrix version of the Chen-Lee-Liu-type DNLS equation 

iQt + Q.. - iQRQ.. = 
iRt - Rxx - iR.QR = 0. 

It is to be noted that there is no restriction on the sizes of Q and R, that is, on ni and 
n2- The system ( |2.5|) was shown to possess at least one higher symmetry |]16|, [l^. The 
system is now proved to be completely integrable in the sense that it has a Lax pair 
and, as a result, an infinite number of conservation laws. 



Next, we shall prove the complete integrability of other systems in [|1^ in the same 
sense. For this purpose, we introduce a transformation of dependent variables: 

Q = F~^QG-^ R = GRF (2.6) 

or equivalently 

Q = FQG R = G-^RF-\ 

Here F and G are invertible matrices, which in general depend on Q and R (or Q and 
R). Then, time evolution equations for Q and R, (|2.5| ), are cast into those for Q and R: 



iQt + Q.X - iQRQcc - iFtF-'Q + iQ{G-%G - {F,F-'),Q + Q{{G-').G}, 
-2F,F-iQ,. + 2Q,{G-%G - 2F^F-^Q{G-%G + {F,,F-^fQ 
+Q{{G-^)^GY + iQRF^F-'Q - iQRQ{G-\G = O 



(2.7a) 



(2.7b) 



iRt - R., - iR^QR - i{G-%GR + iRFtF-^ + {{G-%G},,R - R{F,F-'] 
+2{G-'),GR, - 2R,F,F-' + 2{G-'),GRF,F-' - {{G-'),GyR 
-R{F^F-y + i{G'^)^GRQR - iRF^F-^QR = O. 

A sufficient condition for ( |2.7| ) to be local and closed equations is that F^F'^, FtF~^, 
{G^^)xG and {G~^)tG are expressed locally in closed forms in terms of Q and R, i.e. 
they do not include terms with integrals, infinite sums, etc. We impose this condition 
on F and G in what follows. A closed expression of the Lax pair for the transformed 
system is given by performing the gauge transformation 



^ = g^ g 



G 



Due to this transformation, the linear problem and the Lax pair for ( |2.5|) are changed 
into those for ( |2.7|) : 



U' = g-'Ug-g-'g., 



-h 



+ c 



Q 



R 



-iF^F- 



\RQ - i{G-^).,G 



(2.8) 
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-2/1 



2/0 



2Q 



2R 



-QR 



RQ 



+ C 



12 



21 



+ i 



Vo 



22 



(2.9) 



Here 



^12 = iQ. + -QRQ - iF,F-'Q + iQ{G-%G 

V21 = -ii?^ + \rQR + i{G-^)^GR - iRF^F-^ 

V22 = '^\{RQx - RxQ) + \rQRQ - iRF^F-^Q + i^RQ{G-')^G 

+i^{G-%GRQ-i{G-')tG. 

The above transformation is a powerful tool; it yields new integrable systems of the 
DNLS type by appropriate choices of F and G. To confirm this, we list six illustrative 
examples (a)-(f) with the definition of F and G, the evolution equations for Q and R 
and the transformed Lax matrix U': 



(a) 



F^h 

{G-% = -i^G-'RQ = -i^RQG-' 



{G-')t = G-'{^{RQ, - R,Q) - i^RQRQ} 

= {liRQx - RxQ) - i\RQRQ}G-^ 

iQt + Qxx - KQRQ)x = O 
iRt — Rxx — i{RQR)x — O 



u' = ic^ 



-h 



+ C 



Q 



R 



(2.10) 
(2.11) 



(b) 



G-' 



-i^FQR = -i^QRF 



Ft = f{^(Q,R - QR,) - iiqRQR} 
= {liQxR - QRx) + ilQRQR}F 



iQt + Qxx + iQRxQ + ^QRQRQ = O 
iRt — Rxx + iRQxR ~ -RQRQR — O 



(2.12) 
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(d) 



U' = iC^ 



-h 



+ c 



Q 



R 



+ i 



IRQ 



iQt + Q XX 
iRt — Rxx 



{G-')t = ilmx - QxR) + i^QRQR}G-^ 

= {\(QRx - QxR) + ■i\qRQR - i^g'i?'}^-^ 
iQRQ^ + iQ'^R^ + iQ^QR - ^Q^RQR + ^Q^R" + ^QRQ^R 
iR^QR + iQ:cR^ + \QRR:c + ]^QRQR^ - ^Q'^R^ - ^QR^QR 



U' = iC^ 



-/ 



+ c 



Q 



R 



+ i 



O 



\{RQ + QR) 



G 



-1 



F, = iRQF = ^RQF 



Ft = {^(R.Q - RQ.) + i^RQRQ}F 

= {l(RxQ - RQx) + ■i\rQRQ + iWq^}f 



2 



igi?g^ - iR^Q^ - iRQQ^ + ^i?2g3 _ lgi?2g2 



iRt — Rxx ~ iRxQR — iR^Qx ~ iRxRQ — ^R^Q^ + -R^Q^R 



u' = ic^ 



+ c 



Q 



R 



+ i 



\rq 



\RQ 



Fx = -i-FRQ 



i-RQF 



Ft = F{^(RQ,-R,Q)-i^RQRQ} 

= - RxQ) - i\rQRQ + i\r^Q''}f 

^Qt + - iQRQx - iQRxQ - iQxRQ + iRxQ^ + - RQRQ^ 

+^R^Q^ + ^RQ^RQ = O 

iRf — Rxx — iRxQR — \RQxR — iRQRx ~\~ ^R^Qx ~l~ iRxRQ ~\~ R^QRQ 

-^i?'g' - ^RQR^Q = O 
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U' = 



Q 



(f) 



R 



i-QRF 



O 



(2.19) 



F{^(Q,.R-QRJ-ilQRQR} 

{\{Q.R - QR.) + i\qRQR - i\Q^R^]F 



iQ. + Q.. + ■^QR.Q - .Q^R. - .Q.QR - -Q^R^Q - -Q^RQR 

+\q'R^ + ^QRQRQ = O 



1 



(2.20) 



iRf — Rxx + iRQxR — ^QxR^ — iQRRx + ~RQ^ R? + ~QRQR 



^Q^R^ - ^RQRQR = O 



U' = 



c 



Q 



R 



-\QR 



URQ-QR) 



(2.21) 



I.e. 



For all of the six examples (a)-(f), the compatibility conditions for F and G ^ 

{Fx)t = {F,)x {{G-%h = {{G-')t}x 

can be checked by a straightforward calculation with the help of ( |2.5| ). As is clear from 
the construction of the gauge transformations, Q and R can be rectangular matrices for 
(a) and (b), while Q and R must be square matrices for (c)-(f). 

Comparing the above results with those by Olver and Sokolov, we find that (|2.5| ) 
and (10) in [0, ( CT) and (7) in (|1|) and (12) in fT^, (|2l^) and (14) in 0, 
( CT ) and (16) in 0, (|2l8|) and (15) in 0, (^) and (17) in |T3] are respectively 
identical up to scalings of variables. The system (|2.5| ) was surveyed in a previous paper 
from the viewpoint of the inverse scattering method ||^. The system (|2.10 ) is a well- 
known matrix generalization of the DNLS equation of Kaup- Newell type . For systems 



( CT ), (CT ), (iH) and ( g:20|) , or (12), (14)-(17) in [T^], we have obtained the 
Lax representations by virtue of a non-commutative version of gauge transformations 
for the first time in this paper. 

The systems ( |2.5| ), ( p.l8| ) and ( p.20| ) are interpreted as matrix generalizations of 
( [L.3| ) with a = 1, 6 = 0. The systems ( p.lOj ) and ( p.l6| ) reduce to ( p,.3|) with a = 1, 
5 = —1/4 in the commutative case. The systems ( |2.12| ) and ( p.l4| ) correspond to ( |1.3| ) 
with q; = 1, 5=1/4 when Q and R are scalar- valued. The symmetry approach shows 
that matrix generalizations of ( |1.3| ) are essentially exhausted by the above statement up 
to scalings and the transposition 



||T7[| . It is remarkable that the system (|1.3|) has matrix 
generalizations only for some choices of a and 6. This feature of matrix generalizations 
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may be explained from our approach in the following way. In the case of scalar variables, 
equation (|1.3| ) is generated by the gauge transformation 

q = —;=(] expl — 2i— / rqda;' > 
^/a ^ a J J 

r = ^^rexp<^ 2i— / qrdx'? 
y/a ^ a J J 

for the DNLS equation by Chen, Lee and Liu |]15| 

iq* + ^xx - iqrq^ = 
ir* - r^x - irqr^, = 0. 

Why cannot we generalize this transformation to the matrix case for the continuous 
parameters a and 6 7 As an illustrative example, we choose F and to satisfy 

= i7RQF lim F = / G'^ = I 
where 7 is a scalar parameter. The explicit form of F is given by 
F = S exp {i7 / RQdx'} 

J Xq 

00 

=/+5:(i7)M„. 

n=l 

Here S is the path ordering operator and 

l-X rxi j-Xn-l 

An= dxi dx2--- dx„R(xi)Q(xi) ■ ■ ■ R(x„)Q(x„). 

J Xq J Xq J Xq 

We can calculate the time derivative of F, Ft, with the help of ( |2.5|) . However, it is 
observed that FtF~^ cannot be expressed in a closed form in general, i.e. it includes 
infinitely multiple integrals. An exception is the case of 7 = 1/2, which gives 

Ft = {^{RxQ - RQx) + i^RQRQ}F 

and, as a result, example (d). The example explains why the transformation ( p.6D is 
effective for finite choices of F and G in the matrix case. 

It is a well-known fact that the spatial part of the Lax formulation for the Chen- 
Lee-Liu-type DNLS equation is common with that for the massive Thirring model, i.e. 
they belong to the same hierarchy. Thus, employing an appropriate time part in the 
Lax formulations in correspondence with the Lax matrices U (or U') in this section, we 
can obtain new matrix generalizations of the massive Thirring model (see appendix). 

3. 'C-integrable' systems 

In the previous section, we have verified that all but two systems proposed by Olver and 
Sokolov are 'S-integrable', i.e. they have Lax representations and can be linearized by 
the inverse scattering method. In this section, we show that two systems left for further 
analysis are 'C-integrable', i.e. they can be linearized by a certain type of transformation 
of dependent variables which resembles that in the previous section. 
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The two systems, (11) and (13) in [|T^, are matrix generalizations of ( |1.2| ) with 



a = 0, (3 = 1. We briefly summarize a solution of ( p,.2| ) for investigating (11) and (13). 
The pair of equations (|1.2| ) is rewritten in a linearized form as 

(pe ■'-o )^ - (pe )^^ = 

under the boundary conditions: lim^^a-p p = lima^.^j.Q s = 0. In terms of functions y{x, t) 
and z{x,t), which satisfy a pair of heat equations 

Vt -yxx = zt + z^^ = 

and the boundary conditions, lim^-^a-Q y = lima-^^^o z = 0, the general solution of ( |1.2D is 
given by 



p = ye "-^"=0 ^^'^^ = yh + [a — (3) f zydx'\ 

^ Jxo ' 

,^1>^'^^' = Jl + («-/5) ryzdx'}"^- 

^ J xn ^ 



s = ze 



ioiaj^p (cf. |2| for a = -p) and 

—a r spdx' —a f zydx' 

p = ye -^"=0 = ye •'^o 

a f psdx' a f yzdx' 

s = ze •'''0 = ze •'^o 

for a = p. 

We proceed to solve (11) and (13) in |^ by generalizing the above method to the 
matrix case. We write (11) in |T7| : 

P = P 

* (3.1) 
Sf = — Sxx SSPjjS -|- 2SPS2:. 

Here P is an ni x n2 matrix and S is an n2 x ni matrix. The boundary conditions 
lim P = O lim S = O 

X^XQ X—>XO 

are assumed for convenience. In terms of A defined by 

A, = A{-SP) At = A{SxP - SP,. - SPSP) (3.2) 
the time evolution equation for S in ( |3.1|) is rewritten in a linearized form as 

{AS)t + {ASU = O. 

Here the consistency condition, A^t = Af^, is checked by a direct calculation by use of 
( ^.1| ). We introduce an ni x n2 matrix Y{x,t) and an 77,2 x ni matrix Z{x,t) which 
satisfy a pair of matrix heat equations 

Yt -Y,, = Zt + Z.,x = O (3.3) 

and the boundary conditions 

lim y = O hm Z = 0. 

X^Xq X^Xq 
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The general solution of ( p.3| ) is obtained by means of the Fourier transformation. Thus, 
if we set 

p = Y S = A-^Z 

this gives the general solution of (|3.1|). Due to the relation A^^ = —ZY, we obtain 

A = l2- r Z{x',t)Y{x',t)dx'. 

Here we have assumed the boundary condition, lim^-^^g A = I2, with I2 being the n2 x ^2 
identity matrix. In conclusion, an explicit form of the general solution of is given 
by 

P(x, t) = Y{x, t) S(x, t) = {I2 - r Z{x', t)Y{x', t)dx'y^Z{x, t). (3.4) 

Finally, we shall derive the general solution of the only system left to solve, (13) in 
|[T7| . For this purpose, we set ni = n2 and perform a change of the dependent variables: 

P = AP S = SA-\ 
Noting the fact that ( p.2| ) is rewritten in terms of the new variables P and S as 

A^ = A{-SP) At = A{S^P - SP^ - SPSP + 2S^P^) 
we obtain the evolution equations for P and S: 

= p^^ - 25^p2 _ 2SPP^ + 2SPSP^ - 2S'^P^ 
St = -S^^ - 2S'^P^ - 2S^SP + 2SPS^ + 2SP^S (3.5) 

+2SPS^P + 2S^P'^ - 2S^PSP - 2S^P^S. 
This is nothing but the system (13) in up to scalings of variables. Thus, by virtue 



of the derivation in the above, the general solution of ( |3.5| ), which is an alternative 
expression of (13), is obtained straightforwardly: 

P(x, t) = A-^P = {l- r Z{x', t)Y{x', t)dx'YW{x, t) 

X " X (3-6) 

S{x, t) = SA = {l ~ r Z{x', t)Y{x', t)dx'y^Z{x, t){l - r Z{x', t)Y{x', t)dx'}. 

^ Jxo ' ^ Jxo ' 

It should be noted that all of the matrices in the above expression are square matrices. 
4. Concluding remarks 

In this paper, we have studied matrix-valued systems of the derivative nonlinear 
Schrodinger (DNLS) type. Applying a kind of gauge transformation to a matrix version 
of the DNLS equation of the Chen-Lee-Liu type with a Lax pair P], we have derived Lax 
representations for all but two systems proposed in Hence, these systems can be 



linearized through the inverse scattering formulation and proved to be 'S-integrable' in 
the terminology of Calogero. As has been clarified in section 2, these systems each are 
connected with the others through transformations of the dependent variables. However, 
it is noteworthy that the transformations cannot be written in a closed form in terms of 
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the matrix- valued dependent variables. More explicitly, if F (or G^^) is not the identity, 
we may not express F (or G~^) for the examples (a)-(f) in section 2 without using the 
path ordering operator, infinitely multiple integrals, etc. 

For the two systems in |]T^ left to prove their complete integrability, we have 
shown that both of them are linearizable and connected with each other by a change 
of dependent variables. The transformations which linearize the two systems can be 
explicitly written in a closed form in terms of the variables Y and Z. This enables 
us to obtain the general solutions of the two systems, which directly proves their 'C- 
integrability'. 

To conclude, we have proved for the first time that all the matrix-valued 
systems proposed in can be integrated by the inverse scattering method or the 
transformations of the dependent variables. The dependent variables of the systems 
take values in either square matrices or, more generally, rectangular matrices of arbitrary 
size. However, it is noteworthy that not all of integrable multi-component systems can 
be expressed in terms of matrix variables of arbitrary size. For instance, we consider 
the following Lax pair. 



u = iC 



-h 



-2/1 



+ c 



Q 



R 



2/0 



+ c 



+ 1 



2Q 



2R 

iQ^ + QRQ 



\RQ 
-QR 



(4.1) 



RQ 



+1 

Substituting 



— iRx + RQR 
^{QxR — QRx) + \QRQR 



-,{RQ.-R.Q) + \RQRQ 



(4.2) 



(|1D and 

iQt + Qxx 
iRt — Rxx 



mto ( 
iQRxQ 
iRQxR 



we obtain the evolution equations for Q and R, 
2iQRQx = O 



and the commutation relation, 
[QR, QxR — QRx 



2iRxQR = O 



O. 



(4.3) 



(4.4) 



Setting 



q 



Q = (^1,92, ' 

^ = *(ri,r2, 

which automatically satisfies the constraint ( |4.4| ), we obtain a system of the Kaup- 
Newell-type DNLS equations 0, 

i9t + Qxx - i < q, Vx > q - 2i < q, r > q^. = 

(4.5) 

irt - Vxx - 1< r, q^ > r - 2i< r, q > Vx = 0. 

Here < , > denotes the inner product between vectors. The vector Kaup- Newell system 
( ^75|) is 'S-integrable' because it possesses the Lax pair (|4.1| ) and ( |4.2|) with changing 
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Q ^ q, R ^r. However, the matrix-valued extension of ( |1.3D with a = 1, S = —1/4, 
([4 .31), is not 'S-integrable' in general for Q and R of arbitrary size. 



Acknowledgments 

One of the authors (TT) appreciates a Research Fellowship of the Japan Society for the 
Promotion of Science for Young Scientists. 

Appendix. Massive Thirring model 



In this appendix, we show a list of matrix generalizations of the massive Thirring model, 
which are respectively a member of the DNLS-type hierarchies studied in section 0. 
For this purpose, we consider Lax pairs with the following dependence on the spectral 
parameter (, 



u = ie 



-h 



+ c 



Q 



R 



+ i 



U 



11 



22 



(A.i; 



V 



m 



" -h 


m 




+ i 


" V^ii 


^2 . 




s 


^22 . 



(A.2) 



Here Uu, Vn and U22, V22 are C-independent square matrices whose size are respectively 
rii X ni and ^2 x 'n,2. Q and V are ni x 7^2 matrices. R and S are n2 x ni matrices, m is 
a nonzero constant. We have derived new Lax pairs for several matrix generalizations 
of the DNLS-type equation (|1.3|) in section ^ To obtain matrix versions of the massive 
Thirring model, we have only to change the time part of the Lax pairs as above. In 
this formulation, the new pair of potentials V and S appears. In correspondence with 
the choices of Ujj (j = 1,2), we can determine Vjj (j = 1,2) so that the compatibility 



condition (|2.2|) yields a consistent set of evolution equations. 

We can obtain four matrix generalizations of the massive Thirring model by the 
above-mentioned method. The result of the choices of Ujj, Vjj {j = 1, 2) and the 
corresponding evolution equations is listed as follows. 

(a) Un = U22 = IrQ Vn- 



V22 = o 



Qt - imV + i-VSQ = O 
Rt + imS - i-RVS = O 



V^ - imQ + i-VRQ = O 
S^ + imR - i^RQS = O. 



Un = f/22 = O Vu 



--VS 
2 



22 



-sv 

2 



(A.3) 



(A.5) 



Complete integmbility of derivative nonlinear Schrddinger-type equations 13 
Qt - imV + \^{QSV + VSQ) = O 

Rt + imS - {SVR + RVS) = O 4) 

Vx — imQ = O 
Sx + imR = O. 

(c) Un = U22 = \{RQ + QR) Vn = -lvS V22 = -\VS 

Qt - imV + \^{VSQ - QVS) = O 
Rt + imS - i^{RVS - VSR) = O 
Vx - imQ + i]^V{RQ + QR) = 
Sx + imR - iUrQ + QR)S = O. 

(0) [/n = -^i?Q V2i = Vrx = h,SV-VS) V22 = \SV 

Qt - imV + iiQSV + VSQ - SVQ) = O 
Rt + imS - i^iSVR + RVS - RSV) = O 
Vx - imQ + i^RQV = O 
Sx + imR - i]^SRQ = O. 

In view of the Lax matrix U the cases (a)-(l)) respectively correspond to ( p.5| ), (a), (c) 
and (e) in section 0. Other choices of U obtained in section ^ lead to systems which 
coincide with one of (a)-(l)) up to the exchange of t and x, etc. All of the systems (a)-(c)) 
are shown to be connected with the others through the gauge transformations utilized 
in section |[ 

The system (|A.3| ) has been obtained by the authors in [25|. The others seem to be 
new integrable systems. For the existence of products of the matrices in the evolution 
equations, Q, R, V and S must be square matrices of the same size in (c) and (5). 



(A.6) 
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